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Abstract 

We discuss an experiment conducted by Nesvizhevsky et al. As it is the 
first experiment claimed to have observed gravitational quantum states, it 
is imperative to investigate all alternative explanations of the result. In 
a student project course in applied quantum mechanics, we consider the 
possibility of quantummechanical effects arising from the geometry of the 
experimental setup, due to the "cavity" formed. We try to reproduce the 
experimental result using geometrical arguments only. Due to the influence 
of several unknown parameters our result is still inconclusive. 



davole-8@student.luth.se 
2 chrt ur-8@student.luth.se 
3 hanwik- 7@student.luth.se 
4 hansson@mt. luth.se 



1 



Contents 



1 Introduction 



2 Background 



2.1 The experiment 3 



2.2 Reasons for further investigations 5 



3 Alternative explanations 



.1 Particle in a box 



6 

6 

}.2 Energy eigenvalues! 6 

}.3 Neutron countsj 7 

i.4 Combination! 8 



4 Conclusions and further investigations 



A Simulation of a potential 



8 
8 



2 



1 Introduction 



A wellknown property of quantum mechanics is the quantisation of the energy 
levels of a particle trapped in a potential well. For instance, the electromagnetic 
and the strong nuclear force create different kinds of potential wells and are 
responsible for many observed phenomena in nature, such as the structure of 
atoms and nuclei. This suggests that a splitting of the energy levels should 
also be observed for particles in the Earth's gravitational field. But, since the 
gravitational field is much weaker, the effect should be subtle and hard to detect. 

In a letter to Nature [Nev02|, Nesvizhevsky et al. claim to have observed 
such quantum effects of gravity acting on ultracold neutrons (UCNs). They 
conducted an experiment in which the UCNs were allowed to flow through a 
cavity with a reflecting surface below and an absorber above. By measuring the 
number of neutrons exiting the experimental setup, they claim to have observed 
discrete energy levels. However, in their argument they seemingly disregard 
the modification due to the absorber, stating that it is "sufficiently perfect". 
Further they argue that the discretisation is related to the sudden increase of 
neutrons coming through at distinct widths between the reflecting surface and 
absorber. However, since the UCNs are restricted by both the reflecting surface 
and the absorber, also the geometrical effects should be considered. The results 
might even be explained by geometrical arguments only. 

The aim with this report is to show that by only using geometrical argu- 
ments, the effects observed by the scientist at Grenoble can be explained. 



2 Background 

Here we give a brief review of the experiment reported in [Nev02]. A similar 
experiment was first suggested by V.I. Luschikov and A.I. Frank in 1978 [Fra78|. 
We also discuss inconsistencies of their theoretical analysis and their results, 
giving us some ideas of how to approach an alternative explanation. 



2.1 The experiment 

The experiment was performed with UCNs flowing between a reflecting surface 
below, and an absorber above. The absorber and the "mirror" create a slit 
through which the neutrons pass, eventually reaching a detector at the other 
end of the experimental setup (Fig. [I]). UCNs are essential to the experiment 
since they offer many advantages. For instance, they have an energy of about 
10~ 7 eV, corresponding to a wavelength of ~500 A or a velocity of ~10 m/s, 
allowing them to undergo total reflection at all angles against a number of ma- 
terials. The low energy also allows for high resolution, and since neutrons have 
a lifetime of the order of 900 s, it is possible to store them for periods of 100 s 
or more. This makes UCNs available and suitable for research and experiments 
in fundamental physics. The leading research on UCNs is conducted at LAN- 
SCE, the Los Alamos Neutron Science Center, and at ILL, the Laue Langevin 
Institute, the latter holds the current "worldrecord" density of 41 UCNs/cm 3 . 
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Figure 1: A schematic layout of the experiment [NevOCj. A horisontal beam of 
neutrons is collimated and enters the gap between the absorber and the mirror. 
Neutrons of certain energies eventually reach the detector on the right, while 
others are absorbed. 



Nesvizhevsky et al. argue that when the neutrons are trapped in a potential 
formed by the mirror (an impenetrable "floor") and the Earth's gravity there 
will be a discrete set of possible energy levels (Fig. ^). The four lowest energy 
eigenvalues are Ei = 1.41 peV, E 2 = 2.46 peV, E 3 = 3.32 peV and £4 = 



4.08 peV. For a theoretical treatment of this potential, see [ Flu99 |. The lowest 




Figure 2: The energy versus the absorber height, with the four lowest energy 
eigenvalues indicated. The probability functions are plotted schematically. Our 



calculated eigenvalues are in accordance with the theoretical values in [Nev02|. 



eigenvalue E\ corresponds to a classical height, E\ = mgz, of about 15 jum. This 
leads the group to predict that when the slitopening is less than this height no 
neutron transmission will occur. They argue that if the quantummechanical 
wavefunction has a spatial extension larger than the opening, it will not fit, 
and the neutrons have no chance of reaching the detector. In the experiment 
the group observed a stepwise increase in the number of detected neutrons as 
they increased the slit height. In particular they observed, as predicted, that 
when the slitopening was less than 15 /jm no neutrons reached the detector, but 
that there occurred a sudden increase after 15 /im and another at about 20 /im. 
Their results are shown in Fig. 0. The curves are from the theoretical analysis. 
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Figure 3: Neutron count versus the absorber height in the experiment of 
| Nev02|1 . The dots with error bars show the experimental data. The solid line is 
the expectation from a "classical" analysis, while the dotted curve corresponds 
to a quantummechanical treatment with only the first state taken into account. 
For the dashed curve several states and level populations are considered. 



The solid line is the fully classical treatment, where the neutron throughput 
does not have a lower threshold and is described by 

Ndassicaliz) ~ Z X * , (1) 

where z is the absorber height. Thus the neutron count increases as z 1 ' 5 , since 
as the height is increased a larger spread of velocities of the neutrons is allowed. 
The dotted line is, according to [ Nev02| l , the neutron count when only the lowest 



eigenstate is considered. This is merely a translation of the classical curve with 
the lowest threshold taken into account, 

~ (z - z ) L5 , (2) 

where zq is the absorber height at which the slit becomes transparent to neu- 
trons. This assumes that the neutrons behave "classically" once allowed to 
penetrate the gap. The dashed line is the curve when level populations of all 
four lowest states were considered. 



2.2 Reasons for further investigations 

A first thing to point out is that, the energy eigenvalues have never been mea- 
sured, i.e., all these results are entirely theoretical. So, the only data available 
to us are the neutron counts, N, at the detector as a function of the absorber 
height z. 

The potential used in the analysis is V{z) = mgz, which is the classical 
gravitational potential of the Earth. 
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The absorber is thought to be perfect, removing the wavefunctions (i.e., the 
neutrons) completely when they extend into it. From a classical point of view, 
neutrons with different energies bounce off the mirror. The different energies 
will result in different heights. This is the classical explanation to why the 
higher states are missing. If one then looks at the true quantummechanical 
picture, the neutrons are now described by standing waves in the gravitational 
potential, and there are no "bouncing" particles. Now the absorber must also be 
described by a potential and the classical argument will not be sufficient. The 
absorber could hence give a geometrical explanation for the quantum states. 

The experimental statistics for larger slit widths is insufficient [Nev02|, so 
we need consider, say, the two first steps, making the task of an alternative 
explanation a bit simpler. 



3 Alternative explanations 

In the first two subsections we try to recreate the energy eigenvalues obtained 
for a particle in the Earth's gravitational field using only geometrical argu- 
ments. However, a more careful study of [Nev02|, as well as contact with the 



experimental group, makes us believe that no energy eigenvalues have actually 
been measured. Instead of recreating the energy eigenvalues we tried to explain 
the jumps in the number of detected neutrons. This is presented in subsection 



3.3. 



3.1 Particle in a box 

We start with a potential consisting of two infinite walls, this should be a fairly 
good approximation. The reflecting surface, the mirror, can be seen as an 
infinite wall but the absorber needs more consideration. 

The problem is easy to solve analytically Flu9Sf| [ Gas96 |, and the energies 
are given by 

where a is the box width. Thus, the first energy eigenvalue of a neutron trapped 
in a box with infinite walls and a width of 15 /jm is E\ = 0.9 peV. The 
first energy eigenvalue of a neutron in the Earth's gravitational field, E\g rav — 



1.41 peV [Nev02|, is in the same range. Therefore the assumption that we could 
approximate the experimental setup with a infinite well is not quite correct, 
but it gives us a hint that it might be possible to find a potential that could 
reproduce all the energy eigenvalues of a neutron in a " gravitational field" with 
only geometrical arguments. In the next section we try to modify the infinite 
well with, a better approximation of the absorber. 

3.2 Energy eigenvalues 

The best results are obtained by using a so-called Wood-Saxon potential, to 
approximate the absorber. However, the only way to get the gravity potential 



6 



results was to allow the Wood-Saxon potential to be very close to the gravity 
potential itself. This of course, was not the result we were hoping for. Conse- 
quently we were forced to extend our discussion further. As mentioned before 
we received more information about the experiment and therefore instead of re- 
producing the energy eigenvalues, we only needed to obtain the correct neutron 
counts and this will be described in the next section. 



3.3 Neutron counts 

In order to reproduce the count rate we had to consider the length (10 cm) 
of the experimental setup. We further assumed that a neutron traversing the 
experimental setup is absorbed as 

N(x,z) =N max (z)e' k ^ x , (4) 

where N max {z) is the neutron density at the entrance x = 0, x is the distance 
along the mirror, and k{z) is an absorber parameter related to how much the 
neutron wavefunction is inside the absorber. The crucial condition that must 
be satisfied is of course, that when the distance z is less than 15 /im no neutrons 
should be able to reach the detector, i.e., N(x, z) must be close to zero. 

In a distance of Ax, an amount of 1 — A(z, Ax) neutrons will be absorbed, 
yielding 



Hence we get 



N(x 2 ,z) = (1- A(z,Ax))N( Xl ,z). (5) 



A(z,Ax) = 1 _ e -*(*)A« j ( 6 ) 



where A(z, Ax) is the area of the probability functions inside the absorber. 
The function k(z) can be obtained from the experimental data with the help of 
Equation (|j), giving 



k(z) = 

±J 



N ou t{x,z) 



(7) 



N m ax (^) 

where L is the length of the cavity. If we consider only the lowest state ipi(z), 
we must recreate N out (z,x) as the dotted curve in Fig. |||] See Appendix |A| for 
a simulation considering only the first state. On the other hand, when using 
the four lowest states, the outcome must be the dashed line in Fig. [|. This 
case yields the probability function 

|Vf = CxIVil 2 + C 2 \^ 2 \ 2 + C 3 |^3| 2 + C 4 |^ 4 | 4 , (8) 

with the normalisation condition 

d + C 2 + C 3 + C 4 = 1. (9) 

Since we do not know the level population parameters Cj, we have the freedom 
to choose them to make our model agree with the experimental data. 

J If the transverse neutron temperature is 20 nK as stated in corresponding to 



~ 1 peV, even the simple infinite box potential can explain the first step. The smallest 
separation (a ~ 15 /im) corresponds to the high energy "tail" of the transverse neutron 
energy. 
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3.4 Combination 

There is also the possibility that the gravity potential indeed effects the quan- 
tum states of the neutrons in this experiment. However we would like to, unlike 
the group conducting the experiment, consider the effects from the absorber as 
well. This only changes the appearance of the potential in Fig. |2|, but not our 
principle model. 

4 Conclusions and further investigations 

The main problem in our attempts to explain the results with a geometrical 
model is that N max (z) is unknown, depending on the spread of energies in the 
neutron beam. Other difficulties are finding realistic potentials describing the 
absorber and that, due to a lack of data, our treatment is static. The time 
independence will exclude important phenomena such as tunnelling effects. 

Hence our results are inconclusive. We have not yet determined a potential 
satisfying our requirements and we await the next report from the experimental 
group, hoping it will contain the information we need. 

We propose the following improvements of the experiment: 

• Rotating the experimental setup by 90° keeping everything else, espe- 
cially the transverse neutron energies, constant. If the same result occurs 
it would indicate that the result is due only to the geometry of the exper- 
imental setup as no gravitational quantum states can form in this case. 

• Increasing the length of the cavity. If the output would decrease it might 
confirm our theory of an absorption per unit length. 

• A measurement of where the neutrons strike the detector. The probability 
distribution should reflect \tp\ 2 (since there is a standing neutron wave, 
the neutron is not falling in a classical sense). 

A Simulation of a potential 

We start by simulating a potential describing the absorber, considering only 
the lowest state. With our algorithm we obtain the first eigenfunction (Fig. ||) 
for z = 15 fim. Then we calculate the normalised probability function and its 
area inside the absorber, 

A(15 fim, Ax) ~ 0.0173 area units. 

By assuming L = 10 cm and N max (z) = 0.3 the function k(lh fim) can be 
determined as 

k(l5 /xm) ~ 0.54991, 

where we have obtained N ou t (10 cm, 15 /mi) from the experimental data. Finally 
we can determine Ax as 

Ax ~ 0.0320259 cm. 
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Figure 4: A simulation of a potential describing the absorber with the output 
of the first eigenfunction. 



Now we compare the theoretical A(z, Ax), obtained from the experimental data, 
with the simulated areas which we get with our algorithm, see Table [l|. The 



Slit width (/im) 


Simulated area (a.u.) 


Theoretical area (a.u.) 


20 


0.0252 


0.01245 


30 


0.0031 


0.08333 



Table 1: A comparison of the areas obtained from the theory versus the simu- 
lated ones. 

simulation and comparison of the areas show the difficulties we had with our 
model. With the assumption that N max (z) = 0.3, our area inside the absorber 
decreases faster than the theoretical one when z is increased. We must find a 
suitable potential describing the absorber so that the two areas are the same 
for all z. Also, the cutoff below z = 15 fim must be satisfied. 
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